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ABSTRACT The inclusion of nanoparticles has the potential to improve the thermal efficiency of the base 
fluid. The field of nanofluid (NF) dynamics has attracted important attention due to its extensive range 
of practical uses like fuel cells, solar energy, medication administration, heat transfer, microfabrication, 
coolant applications, and other related domains. The aim of this study is to scrutinize the impact of Lorentz 
force, thermal energy, joule heating, heat source and injection parameters, and Brownian and thermoporetic 
diffusions on the hybrid nanofluid over the moving wedge. The stability inquiry is reported for the existing 
work in order to confirm the stable solutions that make the study unique. Novelty of the existing work is 
to investigate the hybrid nanofluid flow and its stability. The nanoparticles MoS». and Ag are suspended in 
ethylene glycol and water used as host fluids. The numerical solution is obtained from the dimensionless 
first-order differential equations, which are achieved from the basic flow phenomena through similarity 
alteration variables. The influence of emerging factors on flow phenomena is reported via graphs. The 
positive eigenvalues report stable solutions, while the negative eigenvalues designate unstable solutions. 
It is perceived that due to Lorentz force, the rate of the fluid declines while the temperature inside the flow 
channel enhances. The velocity profile decreases while the temperature and concentration increase with 
increasing quantities of permeable factors. Similarly, the Forchheimer number causes to enhance the flow 
rate and decrease the heat and concentration outlines. The current analysis is validated by the published 
work. 


INDEX TERMS Hybrid nanofluid, dual solutions, radiated heat flux, joule heating, moving wedge surface, 
porosity. 


I. INTRODUCTION 

By addition the nanoparticles, the host fluid’s thermal effi- 
ciency may be raised. Because of its numerous practical 
applications—such as in fuel cells, solar energy, medication 
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delivery, heat exchange, microfabrication, coolant, and other 
areas—nanofluid (NF) flow has received a lot of interest. 
The idea of a NF was initially presented by Choi [1] using 
suspended nanoparticles in the base fluid. It was determined 
that heat transmission in NF is developed than the base fluid. 
The chemical process occurring within the nanofluid along a 
horizontal tube was studied by Umavathi and Chamkha [2]. 
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Analytical solution has been obtained. It was investigated 
that the chemical reaction factor destabilized the channel. 
The stability analysis was analyzed by Ketchate et al. [3] to 
investigate the non-Newtonian (NN) movement in the context 
of thermodynamic and magnetic fields. It was concluded 
that first solutions is sure while the second is unreliable. 
Puneeth et al. [4] examined the influence of bio-convection 
on the movement of a pseudoplastic NF over a revolving 
cone in the free stream numerically. It was reported that ther- 
mophoresis factor declines the concentration but improved 
the temperature distribution. Irfan et al. [5] explored the 
phenomenon of thermal relaxation in complex Carreau fluid 
with a heat source. Hussain et al. [6] explored the exam- 
ination of NF motion influenced by velocity slip above a 
nonlinear stretched Riga surface with changing thickness and 
Joule heating. The mobility of the nanofluid in the cooling 
interaction was examined by Rafati et al. [7]. According 
to this analysis, the temperature significantly decreased, but 
the heat rate significantly increased. Razzaq and Farooq [8] 
explored the analysis of non-similar forced convection per- 
taining to the flow of Oldroyd-B fluid through a vertically 
stretched sheet. Numerical results have been obtained by 
using bvp4c approach. It was reported that temperature 
distribution declined with enhancing Prandtl number. Mag- 
netohydrodynamics finds use in a number of fields, including 
nuclear electricity, electronic power plants, MHD engines, 
and electronic turbines. Razzaq and colleagues [9] conducted 
an investigation on the non-analogous resolution pertaining 
to the magnetized movement of a Maxwell nanofluid over a 
surface that is being exponentially extended. Cui et al. [10] 
investigated the thermal effect in bioconvection with 
chemically reaction over horizontal surface of Oldroyd-B 
fluid. Similalry, Lorentz force effect with entropy across 
the trapezoidal cavity comprising NFs was inspected by 
Mahmoudi et al. [11]. Their outcomes displayed that the 
addition of the NPs decreased the generation of entropy. 
Additionally, when the magnetic force became stronger, mag- 
nitude of entropy generation increased. Selimefendigal and 
Oztop [12] quantitatively inspected the magnetized hybrid 
conduction of NFs. This analysis reported that fluid veloc- 
ity decrease with magnetic field. Oztop et al. [13] reported 
the phenomenon of mixed convection in a magnetohydro- 
dynamic NF flow within wavy walls. The heat transmission 
level may be changed by growing or decreasing the pro- 
portion of NPs volume based on Hartmann and Richardson 
numbers. 

Hybrid nanofluids, or HNFs, are frequently used to 
improve thermal performance and, consequently, heat rate 
transmission. Jan et al. [14] evaluated the behavior of MHD 
Sisko NF flow with viscous dissipations. It was reported that 
thermophoresis factor increases the concentration but Lewis 
number has inverse impact on concentration distribution. 

Riaz et al. [15] reported computational analysis of heat 
transfer characteristics in a non-similar ternary HNF flow 
across a linearly extended sheet. In order to investigate 
HNF containing Ag-CuO, Hayat and Nadeem [16] used 
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water as the fundamental fluid throughout a spinning tube. 
A mathematical investigation of the impact of magnetic 
hydrodynamics on 3D HNF via extended sheeting was con- 
ducted by Devi and Devi [17]. It is originated that the melting 
effect is higher once the Fourier number varies from 0 and 0.5. 

Chamkha et al. [18] used a horizontal tube to evaluate the 
efficiency of NF and HNF in order to assess the heat process. 
It was revealed that heat rate for HNF is higher than NF. 
Amala et al. [19] examined HNF through a moving vertical 
frame under the influence of heat sources and hall prospec- 
tive. Analytical solution was obtained through Laplace 
transformation. For the same values of volume fraction, heat 
rate for HNF is higher than NF. Ghadikolaei [20] evaluated 
the enviroeconomic impact of solar PV cell cooling technique 
on carbon dioxide emissions reduction. Ghadikolaei [21] 
investigated how to increase the efficiency of photovoltaic 
solar panels using cooling technologies. Ghadikolaei and 
Gholinia [22] investigated 3D mixed convectional MHD 
movement of GO-MoS2 HNF in H20-(CH20OH) hybrid host 
liquid towards the impact of H2 bond. The same authors 
described the fantastic influence of H2 using convection-free 
motion of HNF in a permeable surface [23]. Similarly, 
same authors [24] considered the natural conduction move- 
ment of MHD caused by MoS2-Ag nanomaterials embed- 
ded in C3H¢02H20 hybrid host fluid with thermal effect. 
Kumaran et al. [25] scrutinized the hydromagnetic driven 
convection movement of a Carreau NF across a wedge/plate 
with plate stagnation. Sudhagar et al. [26] considered the 
non-Darcy impressions of mixed convection NF across a 
wedge in a transparent medium. Gaffar et al. [27] studied 
mathematical exhibiting and outcomes for mixed convec- 
tion BLFs of NF generated by a non-isothermal wedge. The 
system was solved numerically and it as concluded that Brow- 
nian motion has decreasing effect on concentration field. 

In recent decades, researchers have paid close attention 
to the Falkner-Skan equation (FSE). Analytical and numer- 
ical procedures were used to produce the FSE’s exact, 
almost accurate, and approximate findings. Mass transit is 
significantly impacted by a number of manufacturing pro- 
cesses, such as nuclear power plants, the cooling procedure, 
and polyester film, copper sheeting for cooling reasons, 
medical treatment strength, and aerospace engineering. 
Mohana and Kumar [28] investigated the consequences of 
Darcy-Forchheimer time-dependent 3D CdTe-C/H20 HNF 
flow across a stretched sheet via convective transport of 
heat. Padmaja and Kumar [29] investigated the viscous 
dissipation and chemical response impacts on MHD NF 
flowing across a rotating channel using RK4 method. It was 
reported that Eckert number enhances the heat transfer 
rate. Hartree [30] provided closed-form solutions for FSE 
in a BLF. Yacob et al. [31] studied FSE over a fixed or 
dynamic wedge with heat analysis using NFF. Similarly, 
Alam et al. [32] investigated FSE in a magnetised nanofluid 
in a BLF through a moving wedge. Numerical consequences 
have been obtained and it is concluded that heat transportation 
rate enhances with FS parameter. 
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In view of the above literature review, there are few articles 
that focus on the motion of stable hybrid nanofluids through 
a moving wedge using the Falkner-Skan wedge flow frame- 
work in convective boundary conditions. The incorporation 
of nanoparticles has the potential to improve the thermal 
efficiency of the underlying fluid. The field of NF flow 
has garnered noteworthy attention due to its wide range of 
practical uses, including but not limited to fuel cells, solar 
energy, medication administration, heat transfer, microfabri- 
cation, coolant applications, and other related domains. The 
purpose of this study is to employ a modified Buongiorno 
structure to analyse the Falkner-Skan characteristics compu- 
tationally for a non-homogeneous combination of ethylene 
glycol and water. Additionally, molybdenum disulfide and 
silver (Ag) tiny particles will be added through a movable 
transparent wedge sheet. The flow theory also takes into 
account the impacts of joule heating, Lorentz force, Brownian 
motion, thermophoresis features, and radiated thermal flow. 
The bvp4c method is employed to compute the complicated 
differential equations. For validation, the existing outcomes 
are equated with published outcomes and outstanding agree- 
ment is originated. The outcome of emerging features on heat, 
velocity, and concentration are revealed through graphs and 
tables. Since some factors exhibit multiple solutions, a stabil- 
ity evaluation is utilized for the stable solution. So, the novelty 
of the existing work is to investigate the HNF and its stability. 


Il. MATHEMATICAL ANALYSIS 

Considering 2D time-independent movement over a porous 
movable wedge sheet with velocity uy(x) = Uwx™/(1 — a) 
and free-stream velocity Ue (x) = Vex” where @ is constant. 
It is interesting to note that the inequality Uy, < 0 indicates 
moving wedge and the equality U,, = 0 characterizes a static 
wedge. The flow geometry is portrayed in Figure |. The basic 
flow equations also involve the key factors such as Lorentz 
force due to magnetic field, thermal radiations, heat rate, and 
Joule heating. The total wedge angle is Q = mz where 


— 2m 
p= m+1° 


NY 
y 


Hybrid nanofluid flow 


: Momentum boundary layer 
_.......: Thermal boundary layer 
—— : Concentration boundary layer 


FIGURE 1. Schematic diagram. 
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Where, m = | represents the steady movement across the 
rotating wedge close to the stagnating zone, whereas m = 0 
represents the flow through the moving wedge close to the 
plate. The y— axis is vertical to the wedge and the x—axis 
is along to the movable wedge. The MBNM that configures 
the influence of thermophysical properties is the basis of the 
current HNF model. 

The basic flows of equations are [35], [36], [37], [38], [39], 
[40]: 


du dv 
Dede dy: 1 
ox oy @) 
Ohnl 2 ULhnl 2 2 
Br+ ) (u—ue) —F (u- — uz), (2) 
(= PoniKk ° ( :) 
oT oT 1 


(pce, CD #7 \ aT B2 
Ps Ms (Ds ght ) OhntB” 2 (3) 
(PCp) nt dy Too dy J dy (0Cp) nn 
aC =69aC.)—SCOUiDy &T aC 
ua tv—— = + Dp (4) 
ox dy = Too OY oy 


With 
u(x) = EUy, V(X) = vy, T(x) = Ty, C(x) = Cy aty = 0, 


u(x) > u(x) = ax™, T(x) > Ta, C(x) > Cog at y > 00 
(5) 


In above equations u,v signify the velocities, T (temper- 
ature), and C (concentration), Ty, (wall temperature), T., 
(free-stream temperature), the wall NP volume fraction is rep- 
resented by Cy, while the free-stream volume Np is denoted 
by Coo. Here, p is density, jz is the viscosity, , v denotes kine- 
matic viscosity, and o is the thermal conductivity. Radiated 
heat flux is denoted by g,, D7 reports the thermal diffusion, 
and and Dg is Brownian diffusion coefficients. The inertia 
variable is F=(F = 4) and porosity is given by the 
relation K = Kox!~™., and varying magnetic field strength 


isB = Box 2. Also the partial volume fraction are 

denoted by xmos, and Agyos,. The thermo physical quan- 

tities for the present model are given in Tables! and 2. 
Furthermore, the value of g, in Eq. (3) is given by 


= —4229T (Rosseland imation) [25]. Furth 
Gr = —3RK* Gy (Rosseland approximation) [25]. Further- 


more, when we describe the term T* as a linear function of 
T while disregarding the higher-order terms, become T* = 
4T3,T — 3T4, so the equation becomes: 

1607, af 
— 3n* ay 
Presenting the alteration of variables, the scheme of equa- 
tions (2)-(5) become while Eq. (1) satisfied identically 


[36], [37]. 
ed on WY =( 2av| ) (omg) 
y Lm | ‘fm \m+i) *~° 
T = (Tw — Too) O(N) + Too, C 
= (Cw — Coo) (9) + Coo (6) 


qr = 


Nie 
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TABLE 1. Thermophysical properties MoS, Ag and host fluids are [24]. 


Pro Expressions 
per 

ties 

nmi Hnni = = 


[1— Gos, + ole 


Phni Pan = (A -pi + XMoS2PMoS> + XAgPAg 


(p Cy) 


(p65), = (1 =)(pep), + 


XMosz (0¢p) mos, as Xag (Pep) 4, 


Kant 
14+ 3x(Xmos,Kmos, + XagKag z xk) 
Kn} = (mos, Kmos, + Xagkag + 2xk1) ky 
X(Xmos,Kmos, + Xagkag =e xk )ky 
Snnl 
- 3x(Xmos,7mMos, + XagFag — Xi) _ 
Phnl| = (Xmos, mos, + XagIag + 2x01) 


x(Xmos,%™os, + XAglAg — X01) 


So, Eqs. (2)- (5) are distorted to: 


A3M* A 
yr - (As ma +e ae 


Aj 
A At 1 
2 pergite(® + ) 
1 


i te) Ce) 
=0, (7) 
(Aq + Rp) 6” + AsPr0’ + PrNbé’¢' 
+ 30-Rp0* + PrNto” + 30,Rpd0" 
+ 667Rp00” + 367Rp076" + 363>Rpd76” 
+ 63Rp0°0" + A3PrEcM*f” = 0, (8) 
NtO” +Nb@” + PrNbLef ¢’ = 0. (9) 
With 
f(n) = S.f'(n) = €,¢(n) = 1,¢() = 1, at n = 0, 
f'(n) >, O(n) > 0, b(n) > 0, at n> ov. (10) 


The physical factors in above equations are: 
piu, 
(pcp), (Tw — Too)’ 
tDg (Cy — Coo) 
VI : 


Ec(Eckert number) = 


Nb = 
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TABLE 2. Thermophysical characteristics of MoS,—Ag/C,H,0 — H,0 
[24], [33]. 


Physical 
properties 


a | 


Base 
fluid 
Haale C,H,02 
enum 
Silver disulfid — H,0 
e 
Density (~) 10490 5060 peel 
3630 
397.21 
0.387 


1063.8 


Specific heat 

(Cy) 
Thermal 
conductivity 904.4 
(k) 


Electrical 
conductivity 


1 
; ky 
Le(Lewis numbar) = 
B (Pcp), 
Cp), Vv 
Pr(Prandtl number) = elt 
1 
Fp(Forchhei ber) 209 
p(Forchheimernumber) = —— 
VK 
. o1BS 
M(magnetic parameter) = ——— 
pa 
. aKo 
P(porosity parameter) = —. 
2v 
These parameters are expressed as: 
k, 
A= Ennl Ae Phnil Aes Onnl Ape hal 
i pl orl ky 
c 1 
As = (e¢r)imt 9 — 9, -1.M* = M, 
(pcp) j m+1 
F (pcp) 
ie Pp es (11) 
m+ 1 (Pcp), 
The physical quantities of interest are given below: 
1 
Cpr = Reg Cp = AV 2m 4 2f7(0) (12) 
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1 i 
Nur, = Rey 2Nuy = -/" oe + R63) 9’(0), (13) 


Sh, = Rey?" = -Ap \e $/(0). (14) 
Ill. STABILITY ANALYSIS 

Together with the shooting procedure, the resultant nonlin- 
ear system is elucidated computationally using the bvp4c 
technique. For the purposes of this investigation, stability 
evaluation is used since certain variables show dual solutions. 
Suction and injection lead to dual solutions. The approach 
used for stability reasons is explained in [33] and [34]. 
Following is the unsteady flow features (1)-(6). 


du ov 


= 0, 15 
ox oy ye) 
a a 2 

Bg aa A gO nn 07Uu 
at ax dy dx Phnt dy? 


(a B 4 hal Jo Ue) F Gi u2) , (16) 
Phnl Phnik 


oT Z aT ie oT 1 (« a°T a) 
Uu Vv => 
at ax dy (pep) nt ae dye ay 


; (Cen) np (00 Dr a aT | OinB? 
(Pcp) nt dy Too dy dy (Cp) hn 
ac aC =: 8AC”~—*é—éC«=@Y*-« V-TT? a2C 
+u +yv + Dg (18) 


= D 
at ax dy Too Oy? dy? 
For stability purposes, the subsequent conversions are pre- 


sented to convert the scheme to ODEs such as described via 
[33] and [34]. 


Nie 
Nie 


n | ge wv =f 2av} ) oP 
y 2vy f(y, T) m+ 1 , 
T(x, T) = (Tw — Too) 0, T) + Too, 
uex™— 
C(x, T) = (Cw — Coo) O(N, T) + Coo, T = a (19) 


Eq.(20) in Eqs.(16) —(19),we get 
a°f (= ” 1 rv A 0°f 
an? A\ dn Ay On? 


Az a (FY. (AsM* 1 
- Bier (5) # Ay +m) 


Ad ; arf 
+4 (6+) -G+a0) 5, a =O (20) 


076 a0 a0 a 
5 + AsPr—+P ape oe 
an an on 


(Ag + Rb) =a 
: a0\? 070 
+ 30,Rp6* + PrNt (| —) + 36-Rp9—, 
an an? 
a0 020 a0 
+ 607?Rp (= ) + 307 Rp6? —, + 36;Rpd* (= 
on an? an 
470 oF 
+ @3Rp Daa + A3PrEcM* (= , ) 


uf] 
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es 
—(+art) —— =0, (21) 
0°60 Ps a 
wit 8 raat 0. (22) 
an? an? on OT 
with 
OF 
f(n, T) = ares (n, T) = 6, 
n 


O(n, T) = 1, d(H, T) = Lat yn = 0, 
(0 T) >, O(y, T) > 0, b(n, T) > O, atn > oo. (23) 


By using the relation, the stability explanations are inves- 
tigated for time-independent flow f(7) = fo(n),9(™) = 
40 (n) , andd (7) = ¢o (n) that satisfied the boundary con- 
straints as described in [33] and [34]. 

£(n, t) = fo (n) +e" I (my, 7), 

4 (n,T) = % (n) +e" F (n,T), 

6 (n, tT) = go(n) +e "Gy, T). (24) 
Unknown parameters, y is the smallest eigenvalues. In view 
of (24), equations (20)-(23) become 


A3M* 1 A2 
yl” J’ jy” 
( re + 5) ae 


4 A3M* ES 1 
Al P* 
(Aq + Rp) FE” + AsPrF’ + PrNbF'G’ 


+ 36,Rp0 + PrNtF? + 30,Rp00F" 
+ 665,RD00F” + 306,RD05F" + 30°RD05 FF 


A 
4, (B+ FRI" 


Ap : ad 
re (6+ Ft) +y7-—==0, 


Cha 
(25) 


OF 
+ 02RpO3F" — a + A3PrECMM*J” + yF = 0, 


(26) 
Nto" + Nbb” + PrNbLef ¢' 
dG 
+yG—-—=0. (27) 
OT 


The boundary conditions are 
I(q, 7) =S,(, 1) =e, GQ, 5) = 1, 
F(n, tT) = latn =0, 
J'(y, t) >, F(y, t) > 0, Gn, tT) > 0, atn > 00. (28) 
So, equations (25)-(28) become 
A3M* 1 A2 A2 
yl" _ oe J! J" ee Fé J” 
( A <P) ta age) 
Poy cc acer (6+ Fr) + J=0, (29) 
Ay! PS) A, De ie 
(Aq + Rp) F” + AsPrF' + PrNbF’G’ 
+ 30,Rp0* + PrN tF? + 36,Rp0F" 
+ 666,RpOF” + 30?Rp67F" + 303>Rp0?F* 
+ 63RpO°F" + A3PrEcM*J? + yF =0, (30) 
NtF” + NbG” + PrNbLefG’ + yG = 0. (31) 
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Variable declaration 


MZ 
Initialization 


Solution evaluation 
through the d 


FIGURE 2. Flow diagram of bvp4c. a. Dual solution presentation. 


With 


J(q) =8, Ty) =e, Gq) = 1, 
F(n) = latn =0, 
J'(n) —, F(n) > 0, G(n) > 0, at 7 > 0. (32) 


IV. NUMERICAL SOLUTIONS AND VALIDATION 

The complex differential issues (7)—(9) related to the bound- 
ary limits (10) are computationally explained in Matlab by 
the byp4c technique. A three-phase integration approach, 
as described. Collocation polynomials are used to achieve a 
homogeneous 4th-order resolution in an integration interval, 
similar to solution Cl-continuous. In the second step, the 
time structure is divided into sub-intervals by applying a 
collocation approach to create a mesh. The mathematical 
system’s solver ensures that a solution exists. To determine 
the error rate for each sub-interval, a solver is employed. If the 
necessary precision is not achieved, the process is carried out 
again using a modified mesh. The process chart for the bvp4c 
approach is shown in Figure 2. The present study is validtaed 
through the available outcomes and excellent settlement is 
found as illustrated in Table 1. With a tolerance of 10~?, the 
results are said to be quite accurate and boundary condition 
n — oo is taken as n = 7 and shows good convergence. 


V. RESULTS AND DISCUSSION 

PDEs are used to construct a mathematical framework for 
the non-homogeneous combination of 50% water and 50% 
ethylene glycol, along with the inclusion of molybdenum 
disulfide (MoS2) and silver (Ag) nanomaterials. In order 
to examine the features of Falkner-Skan features, the PDEs 
are transformed into a dimensionless set of equations with 
respect to nonlinear ordinary differential equations and com- 
putationally calculated using bvp4c over a moving permeable 
wedge sheet adopting the modified Buongiorno framework. 
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TABLE 3. Assessments of the present and available results when 
Fp = M=0,Py = 0.73,Rp = M = Ec = 0. 


f"(0) —8'(0) 
m Ref. [37] | Present Ref. [38] | Present 
study study 
0.0 0.46960 | 0.46841 0.42015 | 0.42113 
0.1 0.50461 | 0.50441 0.42578 | 0.42578 
0.2 0.56898 | 0.56863 0.43548 | 0.43525 
0.3 0.65498 | 0.65472 0.44730 | 0.44711 
0.5 0.73200 | 0.73212 0.45693 | 0.45673 
0.8 | 0.80213 | 0.80204 0.46503 | 0.46501 
1.0 0.92765 | 0.92752 0.47814 | 0.47805 


The influence of joule heating, suction/injection, magnetic 
field, thermophoresis characteristics, Brownian motion, and 
radiated heat flux on the flow characteristics are elabo- 
rated through graphs and tables. Together with the shooting 
procedure, the subsequent nonlinear differential system is 
elucidated computationally exhausting the byp4c method. 
As demonstrated in Section IV, the reported results are veri- 
fied with established results. For particular physical variables, 
dual solutions arise because of suctions or injections. As a 
result, the stability exploration is applied, as demonstrated in 
Section III. The bvp4c method is used to numerically solve 
the transformed equations (29)—(32) and fixed the parameters 
values as p = 0.1, Ra = 0.2, Fy = 0.1, Np = 0.3, N; = 0.1, 
and Le = 1.1. These values are fixed in the whole investiga- 
tion except when they are shown in graphs and tables. 
Figure 3 provides a graphical representation of the dual 
solutions. According to Wiedman et al. [34] and Merkin [33], 
the suction parameter has an unstable solution for negative 
eigenvalues while being stable for the smallest positive eigen- 
values. The crucial point is the point of convergence when the 
first and second branches come together. When there are mul- 
tiple branches, the stability examination is significantly used 
for analyzing the stable branch. Figure 3 presents a visual 
depiction of the double solutions. The presence of negative 
eigenvalues indicates an unstable outcome, while the smallest 
positive eigenvalues indicate a stable solution. The critical 
point is that when the first and second branches converge. The 
smallest eigenvalues for different values of magnetic factor 
are also depicted in Table 4 when Nb = 0.3, Nt = 0.4, 
Le = 2, 6, = 1.2, P=2, Fp =2, Rp = 1.5, Ec=0.1. 
Figures 4 and 5 characterize the impression of @ag on 
the Skin Friction (SF) and Nusselt Number (NNb) about M, 
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TABLE 4. The first and second solutions for y regarding several 
magnitudes of M when Nb = 0.3, Nt = 0.4, Le = 2, (w= 1.2, P=2, 
Fp = 2, Rp = 1.5, Ec=0.1. 
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FIGURE 4. Inspiration of gag on f” (0). 


Mcl= 0.5454 
Mc2= 0.5035 
Mc3= 0.1716 


FIGURE 5. Inspiration of bag On —6'(0). 


separately. Dual assessment takes place for ag. The hard 
line intelligences the first solution although the dotted line 
pinpoints the unstable solution. The critical magnitudes for 
gag = 0.001,0,01, and dag = 0.1 are Mp) = 0.5454, 
M.2 = 0.5035, and M,; = 0.1716, correspondingly. It is 
worth mentioning at this point that SF boosts with increasing 
magnitudes of ¢4g in the stable solution whereas diminishing 
in the unstable solution, as shown in Figure 4. Moreover, the 
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M 
first solution Second 
solutions 
1 1.4023 —1.2358 
0.8 1.1062 —1.0881 
0.7 0.8505 —0.8245 
0.5 0.6311 —0.6381 
0.4 0.2025 0.3273 
0.2 0.0320 —0.0803 
0.15 0.0061 —0.0002 
FIGURE 6. 
Ec1l= 0.4995 First solution 
€c2= -0.8920 
€c3= -2.8731 
| 


“4 3 2 a 0 1 


N 
w 


FIGURE 7. Inspiration of S on —6’(0). 


heat level drops for the stable solution however increasing 
for unstable solution. The relation between the free stream 
velocity ¢ and the suction on the SF as well as NNb is shown 
in Figures 6 and 7, correspondingly. The point which connects 
the first and second solutions with respect to the e and S are 
Ec] = 0.4995, €-2 = —0.8920, and €.3 = —2.8731 regarding 
to the suction S = 0.5, 1.0, and 1.5, respectively. It is 
investigated that the SF and NNb are increasing function of 
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Scl= -1.7508 
Sc2= -1.9100 
Sc3= -2.1102 


f"(0) 


FIGURE 9. Inspiration of m on —6’(0). 


é and S, in the stable branch but decreasing in the unstable 
branch. In Figure 7 it is noteworthy to footnote that double 
manifestation arises in the unstable branch. In the domain 
—3 <e < 2, deterioration performance is perceived for the 
NNb whereas in the interval 2 < e < 3 growing actions is re- 
counted. Moreover, the impression of rotational angle m on 
the SF with NNb is exhibited in Figures 8 and 9, respectively 
regarding S and Pr. It is remarkable to annotation here that as 
m is increased, the profiles of f” (0) and —6’(0) are sketched 
to be enhanced in the first branch, which gives excellent 
agreement with the results reported by Ibrahim and Tulu [27]. 
This examination demonstrations a 3.8% escalation in heat 
as m is improved. This examination evidenced that once m is 
boosted, the thermal productivity of the base fluid develops. 
Likewise, it is described that at ¢ = 1.0, the SF approaches 
zero i.e., f” (0) = 0, owing to no frictional drag force on the 
sheet while heated convectively. 

The stability analysis has been discussed in detail by 
sketching the above graphs. In the following graphs 9-17, the 
impact of physical factors such as, wedge surface parame- 
ter M, magnetic factor P, permeable factor, Fp Forchheimer 
number, Rg thermal factor parameter, Lewis number, Eckert 
number, Brownian motion factor, and thermophoresis factor 
on the dimensionless velocity, temperature, and concentration 
profiles are discussed in detail. Graphs 18 and 19, report 
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FIGURE 10. f’(y),6(n), and ¢(n) against M. 


the variation of streamlines and isothermal flow for various 
values of magnetic and radiation factors. Similarly, Table 5 
depicts the considerable effect of the nanoparticles on physi- 
cal parameters of interest. 

Figure 10 addresses the influence of magnetic factor M 
on f’(n),0(y), and ¢ (7) profiles. It is significant to report 
here, that profiles f’(7) and 6(7) are growing function of 
M whereas the concentration profile (7) decreases with 
increasing M. The quantities of the magnetic field impression 
denote the concentration of Lorentz forces. The aforemen- 
tioned force is a resistive force which reduces the velocity 
of HNF. However, in the present study the fluid velocity 
raises with the growing magnitude of M. This conclusion a 
correlation with the outcomes of Ibrahim et al. [27]. It should 
be prominent that when M = 0, results in a TBL. 

Figure 11 reveals the impact of m on f’(y),6(n), and ¢ (7) 
distributions of the HNF. It is interesting to note here that as m 
is increased, the profiles of f’ (7), and 0(7), are sketched to be 
enhanced in the first branch and second branch, which gives 
excellent agreement with the results reported by Ibrahim and 
Tulu [27]. This exploration displays a 3.7% escalation in 
heat as m is improved. The investigations demonstrated that 
while m is boosted, the thermal proficiency of the base fluid 
enriches. Likewise, it is described that ate = 1.0, the SF 
turn into zero ie., f’” (0) = 0, owing to no frictional drag 
force on the wedge surface when heated. The effects of m 
on f’(7),6(y), and (ny) profiles are shown in Figure 11. 
In this case, higher values of m fallouts in a reduction in the 
profiles f’(7) while an escalation in 6(7) and (7) profiles are 
observed. The factor m is contrariwise related to 6 for m< 1. 
As the fluid flow slows, an escalation in the wedge angle 
causes a lessening in velocity. This result is consistent with 
Ibrahim and Tulu [27]. The variation of the variable P (perme- 
ability factor) on the profiles f’(7),6(n), and @ (7) is reported 
in Figure 12. It is perceived that the velocity distribution f’(7) 
decelerates as the quantities of P enhances. Physically, when 
the permeability variable grows, the fluid motion become 
slowly, this decays the momentum BLT significantly. Due to 
slow motion of the fluid over the isothermal wedge sheet the 
HNF heated significantly which increases the TBL. 
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FIGURE 12. Variation of f’(7),0(n), and ¢(n) against P. 
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FIGURE 13. Variation of f’(7),0(), and ¢() against Fp. 


Figure 13 has been drafted to see the inspiration of Fy 
on dimensionless velocity, temperature and concentration 
profiles. Form this study dual branches have been observed. 
The fluid velocity enhances when the quantities of F, is 
increased while opposed trend is detected for the 6(7), and 
¢ () profiles. Comparing the three profiles, it is analyzed 
that variation is higher for the second branch than the first 
branch. This analysis is more visible from the velocity field. 

The Brownian motion factor Np effect on 6(7), and #(7) is 
reported in Figure 14. Because of the random motion of the 
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FIGURE 14. Variation of 0(7), and ¢(7) against Np. 
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FIGURE 15. Variation of (7), and ¢(7) against Nt. 


tiny fluid molecules, the energy generated within the fluid 
which enriches the fluid temperature but the concentration 
decreases. Comparing the profiles of 0(7), and ¢(y), the 
unstable solution is greater than the stable. The variation 
of 6(7) and (7) due to themophoresis factor N; is shown 
in Figure 15. When N; increases, due to thermophoresis 
mechanism the heated molecules are moved away from the 
hottest zone to the smallest zone which causes the inflation of 
the HNF temperature for both stable and unstable solutions. 
While, the concentration outlines decrease. 

Figure 16 portrays the performance of 6(7), and ¢(7) when 
the Le is augmented. It is significant to note that both profiles 
decline as the Le is increased. In physics point of view, when 
Le is increased, the mass diffusivity of the fluid decreases due 
to which thermal diffusivity enhances which decline the tem- 
perature and concentration. It is perceived that the increase in 
the first solution is more significant than the second one. 

Joule heating characteristics can be expressed by the 
product of Ec and magnetic factor. Already, discussed in 
aforementioned graph (see Fig. 10) that as M is increased 
the temperature is improved. In the same way the effect 
of Ec on temperature is same. Similarly, thermal radiation 
Rg impact on the temperature field is directly proportional. 
As the quantities of Rq are enhanced the temperature inside 
the fluid increase and this increasing is higher for the second 
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FIGURE 17. Variation of (7), and (7) against Ec. 
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FIGURE 18. Variation of 6(7) and ¢(y) against Rg. 


solution. When the radiation factor is improved, more heat is 
generated due to the radiation reaction in HNF which enhance 
the temperature inside the fluid. 

Streamlines are designed for two dissimilar quantities of 
M = 1.0 and M = 1.2 as publicized in Figure 19 (a-b). 
For dissimilar quantities of M streamlines are dropped. The 
isothermal movement for the assumed model is reconnoitered 
in Figure 20 (a-b) by altered magnitudes of radiated param- 
eter. It is publicized that as the magnitude of Rg upturns 
from Rg = 0.2 to Rg = 1.2, the contour plots of heat 
boosts. Table 5 shows the significance of MoS2 and Ag on 
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FIGURE 19. (a). Streamlines for M = 1.0. (b). Streamlines for M = 1.2. 


TABLE 5. Influence of XMS, > XAg ON Cy, Nur, and Sh; when 
x = 0.20, m= } M = 0.1, P =2, Fp =2, Rp = 1.5, Ec = 0.12, 
Nb = 0.3, Nt = 0.4, Le = 2, 6w = 1.2. 


Sh, Er En Lae 
1.24200675 si: -------- ——— ——— 
1.30510430 04.00% 0.560% 0.535% 
1.32423336 9.7664% 1.200% 1.123% 
1.36000460 15.157% 2.700% 2.574% 
1.40301560 20.403% 4.077% 4.011% 

XMoS2 XAag Crr Nu, 

0.4 0.1 3.06202 | 0.11045 
0.3 0.3 4.30085 | 0.11646 
0.2 0.2 4.60022 | 0.12071 
0.1 0.3 5.06307 | 0.21258 
0.0 0.4 5.13507 | 0.13034 


the improvement rate Ey, E,, and E,, of the physical quantities 
of interest Cy, Nur, and Sh, Additionally, it is shown that the 
related enhancement rates Ey, En, and E,, are significantly 
increased when the solid nanoparticles MoS2 are sus- 
pended in a volumetric percentage of solid nanomaterial Ag. 
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FIGURE 20. (a). Isothermal flow for Ry = 0.2. (b). Isothermal flow for 
Rg = 1.2. 
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TABLE 6. Analysis of —0’(0) for Ag, MoSz, and Ag—MoS). 


Nt | Ag MoS2(—8'(0)) | Ag —MoS2 
(—0'(0)) (—6'(0)) 

0.1 1.33433718 | 1.42221003 1.64443005 
0.3 1.33288700 | 1.41845170 1.6396390 

0.5 1.3314662 1.41470446 1.63690668 


Table 6 represents the reduced Nusselt number for the NF and 
HNF using dissimilar values of Thermophoresis parameter. 
It is observed that HNF has higher heat transmission rate as 
related to NF. 


Vi. CONCLUSION 

A mathematical model is established to scrutinize the effect 
of Lorentz force, thermal radiation, joule heating, heat source 
and injection parameters, Brownian and thermoporetic diffu- 
sions on the hybrid nanofluid over the moving wedge. The 
stability exploration is reported for the existing study in order 
to confirm the stable solution which makes the novelty of 
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the study. The nanoparticles MoS» and Ag are suspended in 

ethylene glycol and water using as a host fluid. The numer- 

ical solution is obtained from the dimensionless first order 
differential equations which are gained from the basic flow 
equations through similarity transformations variables. The 
influences of emerging aspects on heat, velocity as well as 
concentration are revealed through graphs and tables. Since 
some factors exhibit multiple solutions, a stability evaluation 
is utilized for the stable solution. So, the novelty of the 
current study is to investigate the hybrid nanofluid flow and 
its stability. 

The key points of the present study are: 

1. It is reported that the point which connect the points 
of gag = 0.001,0,01, and dag = 0.1 are Moy = 
0.5454, Mc2 = 0.5035, and M,; = 0.1716, correspond- 
ingly. 

2. Itis worth mentioning here that SF boosts with rising mag- 
nitudes of @4g in the stable solution however diminishing 
in the unstable solution. 

3. The point which connects the stable and unstable solutions 
with respect to the ¢ and S are ¢,; = 0.4995, e422 = 
—0.8920, and ¢,.3 = —2.8731 regarding to the suction 
quantities S = 0.5, 1.0, and 1.5, respectively. It is investi- 
gated that the LSF and LNN are the escalating function of 
é and S, in the stable branch but decreasing in the unstable 
branch. 

4. It is interesting to note that as m is increased, the pro- 
files of f” (0) and —0’(0) are sketched to be enhanced in 
the stable branch, which contributes excellent agreement 
with the results reported by Waini et al. [34]. This explo- 
ration demonstrations a 3.7% escalation in heat when m is 
boosted. 

5. The distributions 6 (7) and f’ (7) have a similar response 
when the parameter m increases while reverse behaviour 
is observed for #(7). 

6. When the Lewis number is higher, it makes the thermal 
and mass boundary layer regions smaller. 

7. Direct relationship is established between the parameters 
Nt, Nb and dimensionless profiles 0(7), and (7). 

8. When the radiation factor is improved, more heat is gener- 
ated due to the radiation reaction in HNF which enhance 
the temperature inside the fluid. 

9. It is publicized that as the values of Rg upsurges from 
Ra = 0.2 to Rg = 1.2, the isothermal flow also develops. 

Future suggestion: This model relies entirely on a two- 

dimensional flow setup that incorporates traditional thermal 

and mass boundary situations, which slightly limits its appli- 

cability. However, if the exploration is to be conducted in a 

three-dimensional manner using the passive control approach 

and incorporating a convective heating or melting heat trans- 
mission development, it could prove to be convenient for 
engineering applications. 
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